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HIGHLIGHTS 


•  The  DWT  is  useful  for  analyzing  signal  with  non-stationary  and  transient  phenomena. 

•  The  discharging/charging  voltage  signal  (DCVS)  is  considered  as  original  signal. 

•  Daubechies  dB3  wavelet  is  used  in  the  multi-resolution  analysis  (MRA)  of  the  DWT. 

•  DWT-based  MRA  is  used  to  extract  variable  electrochemical  information  from  the  DCVS. 

•  The  proposed  work  shows  the  clearness  for  reliable  state-of-health  (SOH)  diagnosis. 
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This  paper  introduces  an  innovative  approach  to  analyze  electrochemical  characteristics  and  state-of- 
health  (SOH)  diagnosis  of  a  Li-ion  cell  based  on  the  discrete  wavelet  transform  (DWT).  In  this 
approach,  the  DWT  has  been  applied  as  a  powerful  tool  in  the  analysis  of  the  discharging/charging 
voltage  signal  (DCVS)  with  non-stationary  and  transient  phenomena  for  a  Li-ion  cell.  Specifically,  DWT- 
based  multi-resolution  analysis  (MRA)  is  used  for  extracting  information  on  the  electrochemical  char¬ 
acteristics  in  both  time  and  frequency  domain  simultaneously.  Through  using  the  MRA  with  imple¬ 
mentation  of  the  wavelet  decomposition,  the  information  on  the  electrochemical  characteristics  of  a  Li- 
ion  cell  can  be  extracted  from  the  DCVS  over  a  wide  frequency  range.  Wavelet  decomposition  based  on 
the  selection  of  the  order  3  Daubechies  wavelet  (dB3)  and  scale  5  as  the  best  wavelet  function  and  the 
optimal  decomposition  scale  is  implemented.  In  particular,  this  present  approach  develops  these  in¬ 
vestigations  one  step  further  by  showing  low  and  high  frequency  components  (approximation  compo¬ 
nent  An  and  detail  component  Dn,  respectively)  extracted  from  variable  Li-ion  cells  with  different 
electrochemical  characteristics  caused  by  aging  effect.  Experimental  results  show  the  clearness  of  the 
DWT-based  approach  for  the  reliable  diagnosis  of  the  SOH  for  a  Li-ion  cell. 

©  2014  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

As  one  of  the  most  used  batteries  in  consumer  portable 
electronics  in  recent  decades,  Li-ion  cells  have  become  more  and 
more  favorable  choice  for  portable  device,  power  electronics,  and 
renewable  storage  applications  [1,2].  As  a  result,  in  these  days, 
battery  management  system  (BMS)  of  the  rechargeable  Li-ion  cell 
is  one  of  the  major  concerns  for  electric-powered  transportation 
such  as  electric  vehicles  (EVs)  and  hybrid  electric  vehicles  (HEVs) 
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in  order  to  guarantee  the  overall  system  performance.  Therefore, 
for  achievement  of  an  improved  BMS,  accurate  and  reliable 
knowledge  of  a  Li-ion  cell  is  crucial  to  achieve  these  goals  [3-8]. 
Failure  to  achieve  an  improved  BMS,  leading  to  over-discharging 
and  over-charging  conditions,  may  cause  permanent  internal 
damage  [9,10].  However,  due  to  the  complex  chemical  and 
physical  process  of  the  cells,  the  behavior  of  the  cells  show 
nonlinear  characteristic,  and  thus  is  difficult  to  acquire  optimal 
BMS.  As  a  result,  in  order  to  overcome  these  weak  points,  the 
electrochemical  models  [11-16]  and  the  electrical  circuit  models 
(ECMs)  [17-23]  have  been  representatively  studied  and  applied 
Li-ion  cell  for  efficient  analysis  cell’s  behavior.  The  electro¬ 
chemical  models  enable  an  utilization  of  a  set  of  coupled  non- 
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Fig.  1.  Windows  functions  of  two  signal  processing  transform. 
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Fig.  2.  Discrete  wavelet  transform  (DWT)-based  multi-resolution  analysis  (MRA)  of 
the  signal. 
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Fig.  3.  Down- sampling  of  the  DWT. 


linear  differential  equations  for  description  of  the  pertinent 
transport,  thermo-dynamic  effect,  and  kinetic  phenomena 
occurring  in  the  cell.  In  addition  to  the  relationship  between 
microscopic  quantities,  the  distribution  in  the  cell  can  be  easily 
translated  into  measureable  quantities  (cell  current  and  voltage). 
Unfortunately,  the  electrochemical  models  should  accompany 
partial  differential  equations  with  numerous  unknown  parame¬ 
ters,  which  cause  a  huge  memory  requirement  and  a  heavy 
burden  of  computation,  thus  these  models  are  inappropriate  to 
real  BMS  [24]. 

The  electrical  circuit  models  (ECMs)  have  been  especially 
investigated  for  the  purpose  of  BMS  and  energy  management 
control  of  EVs  and  HEVs.  These  models  are  generally  comprised 
of  resistance-capacitance  (RC)  networks  such  as  first-,  second-, 
and  third-order  RC  ladders  in  order  to  capture  I—V  dynamical 
characteristics  of  a  cell.  Through  physical  parameterization  of  the 
ECM,  the  cell’s  behavior  can  be  better  understood.  However, 
from  the  model-based  state-of-charge  (SOC)  estimation  point  of 
view,  it  has  been  observed  that  difference  in  electrochemical 
characteristics  due  to  manufacturing  variability  and  degradation 
with  use,  result  in  wrong  parameter  values  that  causes  erroneous 
SOC  estimation  and  low  BMS  performance.  Although  the  error 
can  be  temporarily  reduced  by  repeated  measurement  of  model 
parameter  values,  such  as  exercise  would  be  very  inefficient  due 


to  time-consuming  testing,  costly,  and  labor  intensive  for 
obtaining  parameter  values  [24-26].  Therefore,  in  order  to 
reduce  these  shortcomings  and  to  enhance  battery  management 
with  a  requirement  for  more  accuracy  and  representation  of  the 
complex  nonlinear  electrochemical  processes,  there  have  been 
some  approaches  that  deals  with  parameters  variability  of 
model-based  SOC  and  state-of-health  (SOH)  under  different 
experimental  condition  [27,28].  The  authors  in  Ref.  [27]  pre¬ 
sented  a  data-driven  SOC  and  power  capability  estimation  of 
lithium-ion  polymer  battery  using  multi-state  joint  estimator 
with  different  degradation  states.  The  authors  in  Ref.  [28  pre¬ 
sented  an  identical  data-driven  multi-scale  extended  Kalman 
filtering  based  parameter  and  state  estimation  against  different 
battery  aging  levels.  Two  references  show  the  clearness  for 
reliable  SOC  and  SOH  information  with  different  experimental 
condition.  Unfortunately,  so  far,  little  definitive  answer  has  been 
given  to  this  question.  In  addition,  due  to  the  high  complexity 
and  long  run-time  performance  that  accompanied  by  enhance 
model  of  sophisticated  method  such  as  above  references,  an 
heavy  embedded  microprocessor  should  be  more  inevitably 
required  to  provide  accurate  results  in  real-time  battery  appli¬ 
cation,  such  as  real-time  SOC  estimation  and  SOH  prediction. 
Besides,  it  is  very  time-consuming  for  design  of  SOC  and  SOH 
estimation  algorithms  that  satisfy  the  optimum  specification.  It 
should  be  well  considered  how  to  evaluate  and  get  a  compro¬ 
mise  in  the  balance  of  complexity  and  accuracy  of  the  BMS. 
Therefore,  in  this  overall  perspective,  without  direct  an  electro¬ 
chemical  or  electrical  modeling  of  a  Li-ion  cell,  a  novel  approach 
as  one  of  the  key  technologies  of  the  BMS  should  be  newly 
discussed  in  order  to  minimize  aforementioned  problems  in  this 
work. 

This  manuscript  aims  to  present  a  new  approach  for  character¬ 
istic  analysis  of  a  Li-ion  cell  based  on  the  discrete  wavelet  transform 
(DWT)  with  focus  on  minimization  of  the  disadvantages  (model- 
based  characteristic  analysis).  The  DWT  is  becoming  a  powerful 
tool  for  analyzing  the  discharging/charging  voltage  signal  (DCVS) 
with  a  non-stationary  and  transients  phenomena  [29-46].  One  of 
its  features  is  multi-resolution  analysis  (MRA)  with  a  vigorous 
function  of  both  time  and  frequency  localization.  Through  DWT- 
based  MRA  requiring  filtering  and  down-sampling,  the  informa¬ 
tion  on  the  electrochemical  characteristic  of  a  Li-ion  cell  can  be 
extracted  from  the  DCVS  over  a  wide  frequency  range.  During  the 
research  specifically  developed  for  this  methodology,  wavelet 
decomposition  based  on  the  selection  of  the  order  3  Daubechies 
wavelet  (dB3)  and  scale  3  5  as  the  best  wavelet  function  and  the 
optimal  decomposition  scale  is  used  to  implement  as  mother 
wavelet.  The  DWT  decomposes  the  DCVS  into  time  and  frequency 
domains  and  focus  on  short  time  intervals  for  high  frequency 
component  (detail;  Dn)  and  on  long  time  intervals  for  low  fre¬ 
quency  component  (approximation;  An).  Specifically,  for  feature 
extraction  of  the  DCVS  through  DWT-based  MRA,  approximation  A5 
and  D 5  components  finally  decomposed  are  considered  in  this 
work.  This  study  particularly  develops  these  works  one  step  further 
by  showing  the  values  of  the  standard  deviation  for  high  and  low 
frequency  components  extracted  from  variable  Li-ion  cells  with 
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Fig.  4.  Frequency  bands  corresponding  to  the  DWT  signal. 
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Fig.  5.  DWT-based  decomposition  and  reconstruction  processes  with  five  levels. 


Fig.  6.  Samsung  SDI  18650  LiCo02  Li-ion  cell  with  a  rated  capacity  (1.3  Ah). 

different  electrochemical  characteristics  caused  by  aging  effect. 
Consequently,  it  can  be  shown  that  the  experimental  results  clearly 
validated  the  DWT-based  approach  for  the  reliable  diagnosis  of  the 
SOH.  All  experimental  results  were  conducted  on  18650  Li-ion  cells 
that  had  a  rated  capacity  of  1.3  Ah  (Samsung  SDI). 


The  remainder  of  this  manuscript  is  organized  as  follows.  This 
manuscript  is  divided  into  five  parts,  including  this  Introduction 
section.  Section  2  shows  the  theoretical  background  on  classical 
wavelet  transform,  which  includes  continuous  wavelet  transform 
(CWT)  and  DWT.  Specifically,  the  basic  introduction  on  the  DWT- 
based  MRA  is  also  provided  in  this  section.  In  Section  3,  the 
experimental  setup  to  acquire  the  DCVS  of  a  Li-ion  cell  is  described. 
An  approach  of  the  DWT-based  characteristic  analysis  of  a  Li-ion 
cell  is  presented  in  Section  4.  Moreover,  SOH  diagnosis  based  on 
the  DWT-based  approach  and  additional  issues  dealing  with  our 
study  for  EVs/HEVs  applications  are  explained  in  an  orderly 
manner.  In  the  final  section,  some  conclusions  and  final  remarks  are 
given. 
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Fig.  7.  Experimental  setup  for  implementation  of  the  proposed  approach. 
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Fig.  8.  Two  plots  for  analyzing  characteristics  of  a  Li-ion  cell  through  DWT.  (a)  Discharging/charging  current  profile  (scaled-down  HEV  current  profile),  (b)  Discharging/charging 
voltage  signal  data  (DCVS). 


2.  Theoretical  concept  of  the  discrete  wavelet  transform 
(DWT) 

2.1.  Wavelet  analysis-mathematical  representation 

For  non-stationary  signal,  the  Fourier  analysis  is  not  effective 
since  it  transforms  the  signal  into  the  frequency  domain  the  time 
information  is  lost.  This  deficiency  of  the  Fourier  analysis  can  be 
removed  to  some  extent  by  analyzing  a  small  section  of  the  signal  at 
a  time,  called  as  windowing.  This  type  of  analysis,  known  as  the 
short-time  Fourier  transform  (STFT)  [47],  however,  has  the  draw¬ 
back  in  that  the  size  of  the  time  window  is  the  same  for  all  fre¬ 
quencies.  On  the  other  hand,  the  wavelet  transform  involves  a 
varied  time-frequency  window  and  can  provide  good  localization 
property  in  both  the  frequency  and  time  domain.  Window  func¬ 
tions  of  two  signal  processing  transform  are  shown  in  Fig.  1. 

Basically,  a  wavelet  is  a  function  i/<t)eL2(R)  with  a  zero  voltage. 

00 

J  mat  =  0  (1) 

—  00 

Then  the  continuous  wavelet  transform  (CWT)  of  a  signal  x(t)  is 
defined  as 

oo 

CWT  (a,  fa)  =  2=  J  x(W*(^)dt  (2) 

—  OO 

where  \j/(t)  is  called  the  mother  wavelet,  the  asterisk  denotes 
complex  conjugate,  while  a  and  b  (a,  beR)  are  scaling  (dilation  and 
translation)  parameters,  respectively.  The  scale  parameter  a  de¬ 
termines  the  oscillatory  frequency  and  the  length  of  the  wavelet, 
and  the  translation  parameter  b  determines  its  shifting  position. 
During  CWT  analysis,  the  wavelet  is  shifted  smoothly  over  the  full 
domain  of  the  analyzed  signal.  It  thus  calculates  the  wavelet  co¬ 
efficients  at  every  scale,  generating  a  huge  amount  of  data.  Due  to 
the  huge  amount  of  data  generated  through  CWT,  training  classifier 
based  on  its  coefficients  at  different  scales  can  often  become 
cumbersome. 

The  DWT  is  derived  from  CWT,  and  its  advantage  is  that  it  does 
not  shift  and  scale  continuously,  and  can  only  be  operated  in 
discrete  steps.  It  also  provides  sufficient  information  and  offers 


high  reduction  in  the  computational  time.  The  DWT  can  be  defined 
as 

DWT(a,  b)  ==  j  (3) 

—  oo  '  ' 

where  a  and  b  are  replaced  by  2J  and  2J/c.  This  defines  a  dyadic- 
orthonormal  wavelet  transform  and  provides  the  basis  for  MRA. 


2.2.  DWT-based  multi-resolution  analysis  (MRA) 

In  multi-resolution  analysis  (MRA)  process,  any  time  series  x(t) 
can  be  completely  decomposed  in  terms  of  approximations,  pro¬ 
vided  by  scaling  function  (pjj<(t)  and  the  details,  provided  by  the 
wavelet  function  ^(t),  and  are  defined  as  the  following: 

<Mf)  =  2-fy(2-'t-k)  (4) 

tlk(t)  =  2-ty(2-k-k)  (5) 

where  j  and  k  are  integers.  The  scaling  function  is  associated  with 
the  low-pass  filters  with  filter  coefficients  h(n),  and  the  wavelet 
function  is  associated  with  the  high-pass  filters  with  filter  co¬ 
efficients  g(n)  such  that  g(L  -  1  -  n)  =  (-1  )nh(n)  (filter  length  L). 
These  filters  are  constructed  from  the  selected  wavelet  function 


(a)  Scaling  function(dB3)  (b)  Wavelet  function(dB3) 


Fig.  9.  Time-domain  waveform  of  dB3  (Daubechies  wavelet  family),  (a)  Scaling  func¬ 
tion.  (b)  Wavelet  function. 
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(a)  (b) 

Fig.  10.  Two  components  through  five  decomposition  levels  of  the  DWT.  (a)  Approximation  components  A\  ~  A5.  (b)  Detail  components  Di  ~  D5. 


Fig.  11.  Two  coefficients  through  five  decomposition  levels  of  the  DWT.  (a)  Approximate  coefficients  a,-,*  (cAi  ~  cA5).  (b)  Detailed  coefficients  dj j(  (cDi  ~  cD5). 
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Fig.  13.  Three  plots  for  original  signal  DCVS  x(t),  approximation  A5,  and  detail  com¬ 
ponents  Di  ~  D5  through  DWT-based  five-level  MRA. 
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Fig.  14.  Comparison  between  original  signal  DCVS  x(t )  and  the  sum  of  six  components 
(reconstruction;  approximation  A5  and  detail  components  D i  ~  D5). 
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Fig.  16.  Approximation  A5  and  detail  D i  ~  D5  components  decomposed  by  the  DWT- 
based  five  level  (original  signal:  discharging/charging  voltage  with  discharging  current 
of  3  A). 


implementation,  the  relationship  of  the  approximation  coefficients 
and  detail  coefficients  between  two  adjacent  levels  are  given  as 


Qj,k  =  <  X(t),  h(n  -  2k)ai- l,n 

n 

dj,k  =<x(t),'l'ik(t)  >  ^g(n  -  2k)aJ_1.n 


(8) 


(9) 


i pj,k(t)  and  its  corresponding  scaling  function  expressed  as 

Eqs.  (6)  and  (7): 


=  V2Y,mm-n) 

n 

=  V2^2g(n)(p(2t-n) 


(6) 


(7) 


with  ]Th(n)  =  ^2  and  J2sin)  =  0.  The  approximations  are  the 
low-frequency  componerils  of  the  time  series  and  the  details  are 
the  high-frequency  components.  The  MRA  leads  to  a  hierarchical 
and  fast  scheme.  This  can  be  completely  by  a  set  of  successive  filter 
banks  as  shown  in  Fig.  2.  Considering  the  filter  bank 


Decomposition  low-pass  filter  Decomposition  high-pass  filter 
1  o  1 


O  ©  © 


O  o 


(M  Reconstruction  low-pass  filter  Reconstruction  high-pass  filter 

1  r 


?  9 


©  o  °  o 

-1 

6  0 


©  <j 

.  .  T 

°  <b  < 

Fig.  15.  Low/High  filters  coefficients  of  Daubechies  wavelet  dB3.  (a)  Decomposition 
process,  (b)  Reconstruction  process. 


where  aj ^  and  dj^  represent  the  approximation  coefficients  and 
detail  coefficients  of  the  signal  at  level  j,  respectively.  Since  x{t)  has 
a  length  of  2N,  then  there  is  a  maximum  of  N  levels  of  decompo¬ 
sition,  or  a  maximum  N  (J  <  N)  of  applications  of  h(n)  and  g(n). 
Therefore,  the  J-level  DWT  representation  of  a  signal  x(t)  can  be 
defined  as  Eq.  (10). 
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Fig.  17.  Two  coefficients  through  five  decomposition  levels  of  the  DWT  (from  Fig.  16). 
(a)  Approximate  coefficients  (cAi  ~  cA5).  (b)  Detailed  coefficients  dj ,/<  (cDi  ~  cD5). 
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Fig.  18.  Three  plots  for  original  voltage,  approximation  component  A5,  and  detail  component  D5  for  direct  current  internal  resistance  (DCIR)  profiles  from  1  A  through  6  A  at  the 
discharge  region. 


*(t)  =  J2  ajM2  ^(2  J2  djk2  ty(2  Jt-k ) 

k  =  0  v  7  j  =  i  k  =  0  v  / 

(10) 

where  J  (J  <  N)  is  the  number  of  levels  of  decomposition.  Each  set  of 
wavelet  transform  coefficients  ajy  and  dj^  for  1  <j<J,  represents  a 
band-pass  filtered  and  down-sampled  version  (Fig.  3)  of  the  orig¬ 
inal  signal  x(t).  The  decomposition  process  can  be  iterated,  with 
successive  approximations  being  decomposed  in  turn,  so  that  a 
signal  x(t)  is  broken  down  into  many  lower  resolution  components. 
This  is  called  the  wavelet  decomposition  tree.  In  practice,  a  suitable 
level  number  n  is  selected  based  on  the  nature  of  the  signal,  or  on 
some  other  suitable  criterion.  Fig.  4  shows  a  qualitative  represen¬ 
tation  of  the  transfer  function  of  the  filters  used  in  the  transform. 
The  resulting  frequency  bands  are  consecutive;  the  limits  for  these 
bands  depend  on  the  level  of  the  signal  and  on  the  sampling  rate. 
For  the  detail  Di,  the  upper  limit  is  half  the  sampling  rate  C/s/2), 
whereas  the  lower  limit  is  /s/4.  For  a  generic  detail  Dn,  the  upper 
limit  of  its  frequency  band  coincides  with  the  lower  limit  of  the 
band  of  the  previous  detail  Dn_i,  and  its  bandwidth  is  half  the 


bandwidth  of  its  previous  detail.  Therefore,  a  generic  detail  Dn 
contains  the  information  concerning  the  signal  components  with 
frequencies  included  within  the  interval  [2_^n+1|fs,  2 ~nfs]  Hz.  On  the 
other  hand,  the  approximation  An  contains  the  low  frequency 
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Fig.  19.  Two  plots  for  obtaining  three  DCVSs  with  different  time  interval,  current 
magnitude,  and  abrupt  changes,  (a)  Discharging/charging  current  profile,  (b)  Dis¬ 
charging/charging  voltage  signal  data  (DCVS). 
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Fig.  20.  Experimental  procedures  of  three  discharging/charging  current  profiles. 


components  of  the  original  signal  x(£)  included  in  the  interval  [0, 
2“(n+i)/s]  Hz  [47].  For  reference,  the  inverse  DWT  (IDWT)  is  to 
reconstruct  the  decomposed  signal  x(t)  through  the  synthesis  filter 
bank  with  up-sampling  operations.  The  DWT-based  decomposition 
and  reconstruction  processes  are  shown  in  Fig.  5. 

3.  Experimental  cell  and  test  workbench 

Experimental  studies  were  conducted  on  Samsung  SDI  18650 
LiCo02  Li-ion  cells  that  had  a  rated  capacity  of  1.3  Ah  (Fig.  6;  1S1P 
unit  cell).  The  operating  voltage  of  this  cell  was  considered  in  the 
2.8-4.2  V  range  and  the  nominal  voltage  of  the  cell  is  3.7  V.  The 
cell  was  charged  at  a  constant  current  of  0.65  A  (0.5C)  until  the 
voltage  reached  4.2  V.  Following  this,  a  constant  voltage  charge 
maintained  the  voltage  charge  maintained  the  voltage  at  4.2  V 
until  the  current  decayed  to  100  mA.  In  case  of  discharging  con¬ 
dition,  a  constant  current  discharge  at  4  A  (3C)  led  the  cell  to 
approach  the  desired  SOC  level  until  fully  discharged  at  2.8  V.  As 
shown  in  Fig.  7,  the  battery  testing  platform  consisted  of  pro¬ 
grammable  stand-alone  instruments  and  a  computer-control  unit 
(C++  MFC-based  virtual  measurement  system  with  a  PC  as  its 
platform).  Stand-alone  instruments  are  independently  operated 
thus  it  support  a  wide  variety  of  functions  in  the  power  and 
measurement  field.  It  comprised  a  power  supply  for  cell  charging 
(Agilent  E3633A),  an  electric  load  for  discharging  (Hewlett- 
Packard  6050A).  The  supported  operating  modes  include  the 
charge  modes  (CC-CV,  CC)  and  the  discharge  modes  (CC,  CP,  CC/CP 
with  pulse).  Predefined  discharging/charging  current  profiles  were 
run  with  the  experimental  test  bench  through  a  C++  MFC-based 
interface  on  a  PC,  which  controls  the  power  supply  and  elec¬ 
trical  load  and  acquires  measurements  for  a  Li-ion  cell.  The 
experimental  results  can  be  precisely  collected  using  measure¬ 
ment  and  recording  system  through  a  data  acquisition  board 
[General  Purpose  Interface  Bus  (GPIB)  communication]  with  a  high 


accuracy.  The  recorded  data  include  discharging/charging  current, 
terminal  voltage,  SOC  via  Ah-counting,  and  temperature.  Collected 
data  is  used  as  inputs  for  MATLAB/Simulink  S-function  simula¬ 
tions.  BMS  controller  based  on  a  self-made  digital  signal  processor 
(DSP)  controller  using  STM32F105VC  for  real-time  control  was 
implemented.  Cell  performance  largely  depends  on  ambient 
temperature,  therefore  all  experimental  tests  were  performed  at 
25  °C  in  an  electronically  controlled  environmental/humidity 
chamber  (Hitachi  U-6652P-CH3).  For  reference,  in  order  to  record 
the  cell  impedance  behavior  at  various  excitation  frequencies, 
electrochemical  impedance  spectroscopy  (EIS)  was  used  (Zahner 
IM6ex). 

4.  Proposed  approach 

4.1.  Feature  extraction  for  characteristic  analysis  of  a  Li-ion  cell 

An  approach  of  the  DWT-based  characteristic  analysis  of  a  Li-ion 
cell  is  newly  introduced.  For  a  scaled-down  discharging/charging 
current  profile  of  an  HEV,  the  experimental  voltage  data  (DCVS) 
measured  for  1.3  Ah  Li-ion  cell  at  25  °C  is  shown  in  Fig.  8.  The  DWT 
is  applied  to  the  DCVS  with  time-frequency  domain,  where  the 
different  characteristics  of  each  discharging/charging  current  signal 
appear  more  clearly.  For  analyzing  and  evaluating  the  DCVS  with  a 
non-stationary  and  transients,  decomposition  process  of  the  MRA  is 
implemented  for  extracting  information  from  these  phenomena. 
Fig.  8(b)  as  original  DCVS  signal  x(£)  is  passed  through  high  and  low 
pass  filters.  The  filters  of  h(n)  and  g(n)  used  to  decompose  the 
original  signal  x(t)  with  down-sampling  and  filtered  coefficient 
vectors.  The  orthonormal  wavelet  filter  of  the  Daubechies  family 
that  suitable  for  a  wide  range  of  applications  is  selected  as  the 
optimum  wavelet  function.  Daubechies  wavelet,  which  belongs  to 
the  family  of  orthogonal  wavelets  filtering  a  DWT,  are  widely  used 
in  solving  a  board  range  of  problems  due  to  its  orthogonal  property, 
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Fig.  21.  Enlarged  and  specific  experimental  results  to  extract  the  features  A5  and  D5  components  of  the  DCVS  according  to  three  current  profiles,  (a)  Discharging/charging  current  profile  1.  (b)  Discharging/charging  current  profile  2.  (c) 
Discharging/charging  current  profile  3. 
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Fig.  22.  Three  plots  for  distinguishing  electrochemical  characteristic  and  state-of-health  (SOH)  of  Li-ion  cells  through  DWT-based  MRA.  (a)  Discharging/charging  current  profile,  (b) 
Approximation  components  A5.  (c)  Detail  components  D5. 


which  is  potent  for  localization  and  classification  of  disturbance 
[48,49].  From  viewpoint  of  implementation,  Daubechies  wavelet 
was  selected  as  a  good  choice  for  the  DCVS  analysis.  In  this 
approach,  the  decomposition  process,  in  which  the  Daubechies 
basis  of  order  three  (dB3),  in  five  decomposition  levels,  has  been 
applied  to  produce  different  frequency  bands  of  the  DCVS.  Through 
the  extracted  features  from  all  these  frequency  components,  it  is 
possible  to  obtain  good  discriminated  results.  In  Fig.  9,  the  wave¬ 
form  and  frequency  responses  of  the  order  3  Daubechies  (dB3) 
scaling  function  and  wavelet  function  are  shown. 

By  down-sampling  (£2)  of  the  MRA,  the  outputs  from  both  filters 
are  decimated  to  obtain  two  components  of  detail  and  approxi¬ 
mation  at  level  1  (A\  and  Di).  Among  them,  the  approximation 
component  A\  is  sent  to  the  second  state  to  repeat  this  procedure 


(decomposition  for  obtaining  A2  and  D2  at  level  2).  Due  to  five 
decomposition  levels,  it  can  be  shown  that  there  are  each  five 
approximation  (A\  ~  A5)  and  detail  components  (Di  ~  D5)  as 
shown  in  Fig.  10.  Additionally,  approximate  ajy  (cAi  ~  cA5)  and 
detailed  coefficients  djtk  (cDi  -  CD5)  which  decompose  to  five  levels 
by  DWT  are  shown  in  Fig.  11. 

For  reference,  the  decomposition  level  n  is  selected  so  that  detail 
Dn  contains  the  frequency  of  the  supply  source  /  then  the 
approximation  with  the  same  level,  An,  only  contains  signal  com¬ 
ponents  with  frequencies  below/.  The  sampling  frequency  used  for 
capturing  the  DCVS  was  fs  =  5000  sample  s-1.  This  sampling  fre¬ 
quency  should  be  considered  to  achieve  good  results  it  is  set  to  an 
integer  multiple  of  it.  With  consideration  of/=  100  FIz,  the  number 
of  decomposition  levels  is  given  by  Eq.  (11). 
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Fig.  23.  Example  for  obtaining  SOH  information  through  standard  deviations  of  approximation  A5  and  detail  D5  components  for  10  arbitrary  Li-ion  cells  that  have  different 
electrochemical  characteristics,  (a)  Standard  deviations  of  10  approximation  components  A5.  (b)  Standard  deviations  of  10  detail  components  D5. 


n  -  Integer 


'log  (fs/fy 

.  l°g(2)  . 


=  5 


(11) 


The  number  of  the  data  window  32,  which  means  that  the 
number  of  resolution  levels  can  be  up  to  five  levels.  As  expressed  in 
Eq.  (11),  the  optimum  number  of  resolution  levels  in  the  analysis 
was  chosen  to  be  five  levels.  The  frequency  band  represented  by 
each  level  depends  on  the  number  of  resolution  levels  in  addition  to 
the  sampling  frequency.  Applying  DWT  to  the  data  window  of  32 
sample  sizes,  with  sampling  frequency  of  5  kHz,  resulted  in  five 
levels  of  resolutions  with  five  approximations  and  details  whose 
frequency  bands  are  listed  in  Table  1. 

Finally,  the  DCVS  as  original  signal  x(t)  is  decomposed  at  the 
expected  level.  After  filtering,  the  two  components  contain  re¬ 
dundancies  and  it  is  valid  to  down-sample  each  of  the  components 
by  a  factor  of  two,  without  losing  any  information.  According  to 
five-level  decomposition,  the  DCVS  can  be  recovered  in  terms  of 
sub-signals  with  different  scale.  Among  all  approximation 
(Ai  ~  A5)  and  detail  components  (Di  ~  D5)  by  five-level  MRA  of  the 
DWT,  with  the  lowest  frequency  band  A5  and  Di  ~  D 5  with  high 
frequency  bands  are  selected  as  shown  in  Fig.  12  and  expressed  in 
Eq.  (12). 


x(t)  —  Aj  +  EDjOVi  —  A'  +  Dj) 
j<! 

At  five  level  x(t)  =  A5  +  +  D2  +  D3  +  D4  +  D5 


(12) 


where  Aj(t)  and  Dj(t)  express  the  approximation  and  the  detailed 
signal  of  the  Jth  level.  Three  plots  for  original  signal  DCVS, 


approximation  component  A5  and  the  sum  of  all  detailed  compo¬ 
nents  Di  ~  D 5  through  DWT-based  MRA  are  shown  in  Fig.  13.  Fig.  14 
is  additionally  shown  for  comparison  between  original  signal  DCVS 
and  the  sum  of  six  components  (Eq.  (12);  approximation  A5  and 
detail  components  Di  ~  D5).  It  can  be  simply  known  that  there  is 
little  difference  between  original  signal  and  approximation 
component  A5.  In  case  of  the  sum  of  detail  components  D 1  ~  D5,  it 
has  positive  or  negative  values  with  zero  basis.  It  can  be  simply 
shown  that  there  is  little  difference  between  original  signal  and 
approximation  component  A5.  For  reference,  as  expressed  in  Eq. 
(12)  and  shown  in  Fig.  14,  by  the  reconstruction  of  the  inverse  DWT 
(IDWT),  it  goes  without  saying  that  there  is  no  difference  between 
original  signal  x(t)  and  reconstructed  signal  of  approximation  A5 
combined  with  detail  components  D 1  ~  D5.  Lower  level  approxi¬ 
mation  and  detail  components  combine  to  create  higher  level 
components  with  up-sampling  and  filtered  coefficient  vectors. 
Fig.  15  and  Table  2  show  the  DWT  filter  construction  showing 
Daubechies  wavelet  dB3  scaling  function  filter  coefficients  corre¬ 
sponding  low-pass  and  high-pass  filters  for  decomposition  and 
reconstruction  processes. 

Various  discharging/charging  current  profiles  are  applied  to  Li- 
ion  cells  in  order  to  extract  features  such  as  current  magnitude, 
time  interval,  and  abrupt  changes  discharging/charging  condition. 
Wavelet  decomposition  based  of  the  order  3  Daubechies  wavelet 
and  scale  5  is  identically  applied.  By  predetermined  scale  5,  for 
feature  extraction  of  the  DCVS  through  DWT-based  MRA,  approx¬ 
imation  A5  and  detail  components  D 5  finally  are  considered  in  this 
approach.  At  first,  in  order  to  analyze  the  DCVS  characteristics  with 


Fig.  24.  Experimental  DCIR  values  for  10  arbitrary  Li-ion  cells  [Qj. 
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Fig.  25.  Schematic  diagram  for  showing  implementation  process  of  the  DWT-based  approach. 


current  magnitude,  the  direct  current  internal  resistance  (DCIR) 
current  profile  [50]  is  simply  applied  to  an  identical  Li-ion  cell.  After 
fully  charging  SOC  1  (100%)  at  a  constant  current  of  0.65  A  (0.5C), 
followed  by  a  rest  period  (1  h),  a  constant  current  discharge  at  4  A 
(3C)  led  the  Li-ion  cell  at  approach  the  desired  SOC  of  50%.  After  an 
identical  period  of  1  h  for  returning  to  the  electrochemical  equi¬ 
librium  condition,  the  DCIR  current  profiles  from  1  A  through  6  A 
(interval  1  A)  at  the  discharge  region  are  applied  during  10  s.  It  is 
made  up  of  single  repetitions  (1  A  through  6  A).  For  analyzing  ob¬ 
tained  6  discharging  voltages,  decomposition  process  of  the  MRA  is 
implemented.  Fig.  16  shows  the  approximation  As  and  all  detailed 
components  Di  —  D5  with  discharging  current  of  3  A.  In  addition, 
approximate  cAs(=ajj<)  and  detailed  coefficients  cDi  ~  cDs{=djtk) 
are  shown  in  Fig.  17.  Three  plots  for  original  voltage,  approximation 
component  A5,  and  detail  component  D5  are  totally  shown  in  Fig.  18. 
Due  to  little  difference  between  original  voltage  and  approximation 
component  A5,  it  can  be  known  that  electrochemical  characteristics 
following  current  magnitude  for  original  voltage  and  approxima¬ 
tion  component  A5  of  the  DWT-based  MRA  are  approximately 
similar.  As  shown  in  Fig.  18,  in  case  of  the  detail  component  D5,  it 


can  be  specifically  seen  that  marked  difference  between  maximum 
value  of  absolute  values  (positive  values)  and  that  of  absolute  value 
(negative  value)  with  zero  basis  is  increased  with  the  magnitude  of 
current  for  identical  discharging  or  charging  time.  Consequently, 
the  approximation  A5  and  detail  components  D5  are  used  to  analyze 
electrochemical  characteristics  with  various  current  magnitudes 
applied  to  a  Li-ion  cell. 

As  aforementioned  before,  for  feature  extraction  of  the  DCVS 
through  DWT-based  MRA,  approximation  A5  and  detail  D5  com¬ 
ponents  finally  decomposed  are  considered  in  this  work.  As  shown 
in  Fig.  19,  three  discharging/charging  current  profiles  including 
different  time  interval,  current  magnitude,  and  abrupt  changes  are 
applied  for  obtaining  each  DCVSs.  Each  flow  chart  for  describing 
experimental  procedures  of  three  discharging/charging  current 
profiles  is  shown  in  Fig.  20.  The  range  of  all  current  profiles  is  -6  A 
through  6  A  with  an  interval  of  1  A.  In  case  of  the  first  current 
profile  of  Fig.  20(a),  constant  discharging  and  charging  currents  are 
sequentially  changed  and  applied  during  1  s  through  10  s  with  an 
interval  of  1  s.  It  is  made  up  of  single  repetitions  (1  A  through  6  A). 
There  is  no  rest  between  previous  discharging  and  next  charging. 
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Table  1 

Frequency  bands  obtained  by  decomposition  in  multi-levels. 


Level 

number 

Wavelet 

decomposition 

Component  type 

Frequency  band  (Hz) 

1 

D, 

Detail 

1250-2500  C/s/22  -  fsl 21)  Hz 

2 

d2 

Detail 

625—1250  (/s/23  ~  /s/22)  Hz 

3 

d3 

Detail 

312.5-625  CL/24  -  /s/23)  Hz 

4 

d4 

Detail 

156.25-312.5  C/s/25  -  /s/24)  Hz 

5 

D5 

Detail 

78.125-156.25  C/s/26  -  /s/25)  Hz 

5 

a5 

Approximation 

0-78.125  (0  ~  /s/26)  Hz 

However,  as  shown  in  Fig.  20(a)  and  (b),  the  second  current  profile 
has  only  rest  time  of  3  s  between  previous  discharging  and  next 
charging.  The  third  current  profile  of  Fig.  20(c)  is  slightly  different 
from  Fig.  20(a)  and  (b).  First  of  all,  DCIR  current  profiles  from  1  A 
through  6  A  (interval  1  A)  at  the  discharge  regions  are  applied 
during  1  s.  After  10  s  for  rest  time,  DCIR  current  profiles  are  iden¬ 
tically  applied  at  the  charge  regions.  It  is  made  up  of  single  repe¬ 
titions  with  duration  for  discharging  or  charging  (1  s  through  10  s). 

Fig.  21  shows  more  enlarged  and  specific  experimental  results  to 
extract  the  features  A5  and  D5  components  of  the  DCVS  according  to 
three  current  profile’s  characteristics  such  as  current  magnitude,  time 
interval,  and  abrupt  changes  discharging/charging  condition.  It  can 
be  known  that  high  frequency  resolution  of  transformed  signal  can 
detect  transient  better  than  the  low  frequency  resolutions.  It  can  be 
seen  that  there  are  a  number  of  sharp  spikes  appearing  immediately 
during  the  period  of  the  transient.  For  long  hours  of  discharging  or 
charging  with  no  change,  there  is  little  difference  between  original 
signal  x(t)  and  approximation  component  A5,  then  the  sum  of  detail 
components  Di  -  D5  is  near  to  zero.  This  effectively  means  that  all 
features  for  an  original  signal  can  be  almost  obtained  from  approxi¬ 
mation  component  A5.  Abrupt  change  of  the  experimental  condition 
of  discharging/charging  current  profile  leads  to  an  increased  sharp 
spike  as  high  frequency  component  D5  (discharging  =>  charging  con¬ 
dition,  charging  =>  discharging  condition).  Namely,  through  DWT- 
based  MRA,  high  frequency  components  Dn  within  a  small  time 
window  can  be  resolved,  and  only  low  frequency  component  An  re¬ 
quires  large  time  windows.  A  high  frequency  component  Dn  com¬ 
pletes  a  cycle  in  a  much  shorter  interval  and  it  is  possible  to  complete 
a  low  frequency  component  a  cycle  in  a  large  interval.  Therefore,  slow 
and  fast  varying  of  the  DCVS  can  be  identified  over  long  and  short 
intervals,  respectively. 

4.2.  State-of-health  (SOH)  diagnosis  through  DWT-based  MRA 

This  present  approach  develops  the  DWT  investigations  one 
step  further  by  showing  approximation  A5  and  detail  D5  compo¬ 
nents  extracted  from  variable  Li-ion  cells  with  different  electro¬ 
chemical  characteristics  caused  by  aging  effect.  For  this,  a  calendar- 
life  test  was  used  to  degrade  the  Li-ion  cell  performance  over  pe¬ 
riods  of  50—150  days  at  60  °C.  When  another  current  profile  shown 
in  Fig.  22(a)  is  applied  to  four  Li-ion  cells,  each  approximation  A5 
and  detail  D 5  components  are  shown  in  Fig.  22(b)  and  (c),  respec¬ 
tively.  The  order  3  Daubechies  wavelet  and  scale  5  is  identically 
applied  for  wavelet  decomposition.  It  can  be  shown  that  voltage 
variance  among  four  cells  are  clearly  related  to  aging  time,  then  the 


magnitude  of  voltage  variance  of  the  cell  degraded  by  calendar-life 
test  during  150  days  is  the  largest  difference  when  compared  to 
that  of  other  cells.  It  is  apparent  that  this  difference  among  the  cells 
can  be  effectively  used  as  the  key  feature  to  distinguish  the  elec¬ 
trochemical  characteristic  and  to  diagnose  SOH  from  the  approxi¬ 
mation  A5  and  detail  D5  components.  Therefore,  through  statistical 
analysis  of  standard  deviation  considering  voltage  variance,  it  is 
possible  to  diagnose  SOH  of  an  arbitrary  Li-ion  cell.  For  reference, 
Fig.  23  shows  the  example  for  obtaining  SOH  information  through 
standard  deviations  of  approximation  A5  and  detail  D5  components 
for  10  arbitrary  Li-ion  cells  that  have  different  electrochemical 
characteristics.  As  shown  in  this  figure,  the  cells  with  the  largest  or 
smallest  values  indicate  the  minimum  or  maximum  SOH  infor¬ 
mation  in  battery  group,  respectively.  Finally,  the  SOH  can  be 
diagnosed  as  follows: 


y^curr  y^aged 

+ 

D^rr_£>aged 

^fresh_^aged 

Dfresh_Daged 

_  soh„5  +  SOHDs 

arbitrary  2  2 

where  A§urr  and  D§urr  are  the  standard  deviation  of  approximation 
As  detail  D5  components  for  an  arbitrary  Li-ion  cell  respectively, 

^feh^fresli)  and  >1aged(Daged) 

are  each  standard  deviations 
approximation  As  detail  D5  components  for  fully  fresh  and  aged  Li- 
ion  cells  in  1.3  Ah  cell  group.  For  reference,  the  fully  fresh  cell  has 
the  largest  SOH  (lowest  direct  current  internal  resistance  (DCIR) 
and  voltage  variance).  On  the  other  hand,  the  fully  aged  cell  has  the 
smallest  SOH  (largest  SOH  and  voltage  variance).  Experimental 
DCIR  values  at  SOC  50%  for  10  Li-ion  cells  are  shown  in  Fig.  24.  Pulse 
power-based  SOH  prediction  of  an  arbitrary  Li-ion  cell  is  generally 
implemented  using  obtained  DCIR  value.  As  shown  in  Fig.  24,  it  has 
been  found  that  the  overall  trend  of  change  between  DCIR  among 
10  Li-ion  cells  is  approximately  similar  with  those  of  approximation 
As  and  detail  D5  components.  Therefore,  it  is  possible  to  obtain  SOH 
information  of  an  arbitrary  Li-ion  cell  using  only  approximation  As 
and  detail  D5  components. 

As  expressed  in  Eq.  (13),  a  final  SOH  of  an  arbitrary  Li-ion  cell  is 
the  average  of  SOHAs  and  SOHDs  which  represent  SOH  information 
for  approximation  As  and  detail  D5  components,  respectively.  Ac¬ 
cording  to  the  values  of  standard  deviation  of  A5  and  D5  compo¬ 
nents,  10  Li-ion  cells  have  each  standard  deviation  of  As  and  D5 
components  within  these  ranges  of  0.1234-0.1593  V  and  0.0676- 
0.0877  V,  respectively.  With  obtained  standard  deviation  values  of 
As  and  D5  components  (0.1404  V  and  0.0769  V)  for  an  arbitrary  Li- 
ion  cell,  three  values  of  SOH^,  SOHDs,  and  diagnosed  SOH  are 
expressed  in  Eq.  (14)  through  (16),  respectively. 


sohAs 


SOHd5 


y^curr  _  y^aged 

0.1404-0.1593 

y^fresh  _  y^ged 

0.1234-0.1593 

D§urr  -  Dl|ged 

0.0769  -  0.0877 

Dfresh  _  Daged 

0.0676-0.0877 

=  0.5264 


=  0.5373 


(14) 


(15) 


Table  2 

DWT  filter  construction  showing  Daubechies  wavelet  dB3  scaling  function  filter  coefficients  corresponding  low-pass  and  high-pass  filters  for  decomposition  and  recon¬ 
struction  processes. 


Low-pass  filter  coefficients 

High-pass  filter  coefficients 

Decomposition 

Reconstruction 

0.0352,  -0.0854,  -0.1350,  0.4599,  0.8069,  0.3327 

0.0332,  0.8069,  0.4599,  -0.1350,  -0.0854,  0.0352 

-0.3327,  0.8069,  -0.4599,  -0.1350,  0.0854,  0.0352 
0.0352,  0.0854,  -0.1350,  -0.4599,  0.8069,  -0.3327 
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S0H  =SOH,5  +  SOHD5  =  0.5264  + 0.5373  =  05318 

arbitrary  2  2 

(16) 

It  is  determined  that  the  range  of  SOH  is  from  0  to  1,  where  1 
means  that  the  Li-ion  cell  is  totally  healthy,  or  it  is  new,  and 
0  means  the  Li-ion  cell  cannot  meet  the  power  demand  of  the 
battery  group.  For  reference,  according  to  total  number  of  experi¬ 
mental  cells,  diagnosed  SOH  range  may  be  changed.  Because,  it  is 
reasonably  certain  the  two  values  of  SOH^5  and  SOHD5  for  fully  fresh 
cell  and  fully  aged  cell  will  be  changed.  For  example,  if  there  are  30 
arbitrary  Li-ion  cells  in  cell  group,  the  cells  with  the  largest  and 
smallest  A^resh(D^resh)  and^5§ed(D5ged)  can  be  used  as  fully  fresh  cell 
and  fully  aged  cell,  respectively. 


4.3.  Additional  issues  considering  DWT-based  approach  for  EVs/ 
HEVs  applications 

According  to  the  proposed  approach,  elaborate  information  on 
electrochemical  analysis  for  a  Li-ion  cell  can  be  provided.  It  is 
possible  to  develop  this  approach  one  step  further  by  showing 
additional  issues  for  EVs/HEVs  applications.  For  electric-powered 
transportation  of  EVs/HEVs,  multiple  cells  are  generally  assem¬ 
bled  together  in  three  configurations  of  series,  parallel,  and  series/ 
parallel  to  satisfy  the  high  voltage/energy  requirements  [51-54]. 
Unfortunately,  because  of  degradation  with  use,  manufacturing 
variability,  and  cell  architecture,  it  has  been  observed  that  indi¬ 
vidual  cells  in  a  real  pack  will  typically  show  some  variation, 
namely  cell-to-cell  variation,  in  performance  depending  on  the 
specific  operating  condition  [55,56].  Inhomogeneity  among  cells 
leads  to  dissimilarity  in  certain  properties  of  the  cells  and  results  in 
performance  variation.  The  non-uniformity  observed  in  cell  volt¬ 
ages  may  limit  the  maximum  output  of  the  battery  pack.  Therefore, 
on  the  pack  level,  the  needs  for  identification  and  equalization  of 
cell  variations  are  importantly  considered  as  critical  issues  for  EVs/ 
HEVs  battery  packs  safety  and  efficiency. 

Various  ways  for  equalization  circuits  have  been  widely  investi¬ 
gated  and  used  to  protect  against  abrupt  failures  to  extend  cell’s  life 
[57—59].  However,  there  are  some  inevitable  limitations  for  over¬ 
coming  cell-to-cell  variations.  Firstly,  it  is  expected  to  be  cost  more 
expensive  and  sizing-increasing  for  achievement  of  the  sophisticated 
equalization  topology.  Secondly,  these  equalizers  cannot  absolutely 
guarantee  SOC  matching  among  the  single  cells,  together  with 
voltage  matching.  Thirdly,  even  though  information  on  voltages  and 
SOCs  of  the  cells  in  the  pack  are  approximately  similar  through 
voltage/SOC  equalizers,  it  is  inevitable  that  voltage/SOC  mismatching 
in  the  battery  pack  certainly  reoccurred  following  the  discharging/ 
charging  current  applied  to  the  battery  pack.  Therefore,  without 
equalization  techniques,  a  novel  approach  should  be  investigated  to 
alleviate  the  problem  of  non-uniformity  caused  by  cell-to-cell  varia¬ 
tion  because  of  electrochemical  characteristics  difference. 

Screening  process  is  to  select  that  have  similar  electrochemical 
characteristics  within  a  cell  group.  Conventional  method  [60] 
implemented  two  kinds  of  screening  processes,  such  as  discharge 
capacity  screening  and  resistance  screening  in  an  orderly  manner. 
The  cells  finally  selected  are  used  to  configure  variable  battery  packs 
with  safety.  Using  these  packs,  Ref.  [60]  describes  an  effort  to  provide 
each  equivalent  circuit  model  for  multi-cell  battery  strings  in  series, 
parallel,  and  series/parallel  connections.  Simplified  model-based 
SOC  estimation  using  extended  Kalman  filter  (EKF)  [60,61]  enables 
us  to  provide  significant  SOC  improvement  with  high  accuracy. 

Through  using  DWT-based  approach,  it  is  more  expected  to 
obtain  elaborate  and  useful  information  on  electrochemical  char¬ 
acteristics  than  conventional  method  of  Ref.  [60].  Above  all, 


discharging/charging  voltage  with  non-stationary  and  transient 
phenomena  can  be  effectively  analyzed.  Besides,  DWT-based  MRA 
is  used  for  extracting  information  on  the  electrochemical  charac¬ 
teristics  in  both  time  and  frequency  domain  simultaneously.  By 
additional  statistical  analyses  including  average,  standard  devia¬ 
tion,  and  maximum/minimum  values  of  the  low  and  high  frequency 
components,  cell’s  similarity  accuracy  of  screening  process  can  be 
increased.  The  BMS  controller  applied  for  real-time  control  such  as 
cell  voltage/current  sensing,  pack  total  voltage  sensing,  and  relay 
control,  is  additionally  used  to  perform  statistical  analysis.  For 
reference,  in  the  proposed  approach,  BMS  controller  based  on  self- 
made  DSP  controller  using  STM32F105VC  is  implemented.  Stable 
configuration  of  the  battery  pack  based  on  screening  process  is 
really  attractive  for  EVs/HEVs  manufactures.  Besides,  outstanding 
results  of  model-based  SOC  estimation  and  SOH  prediction  can  be 
expected  to  provide  en  efficient  EVs/HEVs  operation  with  optimum 
BMS.  In  the  light  of  these  considerations,  DWT-based  approach  is 
absolutely  necessary  and  indispensable  part  of  car  industry.  The 
schematic  diagram  that  depicts  implementation  process  of  the 
proposed  approach  including  additional  issues  for  EVs/HEVs  ap¬ 
plications  is  displayed  in  Fig.  25. 

5.  Conclusion 

Accurate  and  reliable  knowledge  of  the  Li-ion  cell  is  absolutely 
necessary  to  satisfy  for  an  improved  BMS.  Nowadays,  the  electro¬ 
chemical  models  and  ECMs  have  been  representatively  researched 
and  applied  Li-ion  cell.  However,  in  case  of  the  electrochemical 
models,  due  to  difficulty  in  implementing  of  partial  differential 
equations  with  numerous  unknown  parameters,  these  models  are 
inappropriate  to  real  BMS.  RC-network  based  ECMs  enables  us  to 
capture  I—V  dynamical  characteristics  of  a  Li-ion  cell.  Unfortu¬ 
nately,  it  is  inevitable  that  difference  in  electrochemical  charac¬ 
teristics  due  to  manufacturing  variability  and  degradation  with  use, 
result  in  wrong  parameter  values  that  causes  erroneous  SOC  esti¬ 
mation  and  low  BMS  performance.  Thus,  in  order  to  minimize 
aforementioned  problems  of  electrochemical  models  and  ECMs,  a 
novel  approach  should  be  absolutely  discussed. 

This  research  deals  with  an  innovative  approach  to  analyze 
electrochemical  characteristics  and  SOH  diagnosis  of  a  Li-ion  cell 
based  on  the  discrete  wavelet  transform  (DWT).  The  DCVS  with 
non-stationary  and  transient  phenomena  is  applied  as  original 
voltage  x(t)  of  DWT.  Specifically,  in  order  to  extract  information  on 
the  electrochemical  characteristics  in  both  time  and  frequency 
domain  simultaneously  from  the  DCVS,  DWT-based  MRA  is  applied. 
As  a  result,  through  using  the  MRA  with  implementation  of  the 
wavelet  decomposition  (order  3  Daubechies  wavelet  (dB3)  and 
scale  5),  it  is  possible  to  extract  the  information  on  the  electro¬ 
chemical  characteristics  over  a  wide  frequency  range.  This  research 
particularly  develops  these  investigations  one  step  further  by 
showing  low  frequency  component  (approximation  An)  and  high 
frequency  component  (detail  Dn)  extracted  from  variable  Li-ion 
cells  with  different  electrochemical  characteristics  caused  by  ag¬ 
ing  effect.  Experimental  results  show  the  robustness  of  the  DWT- 
based  approach  for  the  reliable  SOH  diagnosis  for  a  Li-ion  cell. 

Nomenclature 

EVs/HEVs  electric  vehicles/hybrid  electric  vehicles 

SOC  state-of-charge 

SOH  state-of-health 

BMS  battery  management  system 

CWT  continuous  wavelet  transform 

DWT  discrete  wavelet  transform 

MRA  multi-resolution  analysis 

DCVS  discharging/charging  voltage  signal 
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An  approximation  component 

Dn  detail  component 

scaling  function  of  the  signal  at  level  j 
i pj,k(t)  wavelet  function  of  the  signal  at  level  j 

ajtk  approximation  coefficients  of  the  signal  at  level  j 
dj^  detail  coefficients  of  the  signal  at  level  j 
h(n)  filter  coefficient  of  the  low-pass  filters 

g(n)  filter  coefficient  of  the  high-pass  filters 

J  decomposition  level 

N  maximum  of  the  decomposition  level 

CC  constant  current 

CV  constant  voltage 

CP  constant  power 

DCIR  direct  current  interval  resistance 
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